The Intermediate Disorder Regime for Directed Polymers in Dimension 1 + 1 



Tom Alberts,^ '[2] Kostya Khanin,^^[^ and Jeremy Quastel^'j^ 

^Department of Mathematics, University of Toronto, Toronto, ON, Canada 
(Dated: March 30, 2010) 

We introduce a new disorder regime for directed polymers with one space and one time dimension 
that is accessed by scaling the inverse temperature parameter f5 with the length of the polymer n. 
We scale /3„ := /3n~" for a > 0. This scaHng sits in between the usual weak disorder (/3 = 0) and 
strong disorder regimes (/3 > 0). The fluctuation exponents ^ for the polymer endpoint and x for the 
free energy depend on a in this regime, with a = corresponding to the usual polymer exponents 
^ — 2/3, X = 1/3 and a > 1/4 corresponding to the simple random walk exponents ^ = 1/2, x = 0. 
For a £ (0,1/4) the exponents interpolate linearly between these two extremes. At a — 1/4 we 
exactly identify the limiting distribution of the free energy and the end point of the polymer. 
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INTRODUCTION 

Directed polymers in disordered media is a model of a 
variety of physical phenomena. It describes vortex lines 
in superconductors [1], domain walls [2], roughness of 
crack interfaces 0, Burgers turbulence [1] and the KPZ 
equation |12| . It provides a simple model in which the 
more difficult predictions of spin glasses can be tested 
and is related [5] to the phase boundary in the 2D Ising 
model. Starting with [2l|5], they have been the subject 
of intensive study O [7| over the last twenty- five years. 
In the setting of the d-dimensional integer lattice, the 
polymer measure is a random probability measure on the 
set of paths of d-dimensional nearest neighbour lattice 
walks. The randomness of the polymer measure is inher- 
ited from a random environment which manifests itself as 
an i.i.d. collection of random variables placed on the sites 
of Z_|_ X Z'^. Given a fixed environment lu : Z+ x Z"* — > M, 
the energy of an rt-step nearest neighbour walk s is 

The polymer measure on such walks is then defined in 
the usual Gibbsian way by 



where /3 > is the inverse temperature, P is the uniform 
measure on n-step walks, and Z!^{/3; x) the point-to-point 
partition function 

Z-(/3;x) = Ep [e-^^"(^)(5o(5(n) = x)]. 

The partition function Z!^{j3) is the sum of the point- 
to-point versions over x. The random environment is 
a probability measure Q on the space of environments 
n = {uj : Z+ X Z"^ ^ E}. Let Q be a product measure 
so that the variables uj{i,z) are independent and identi- 
cally distributed, and assume they have mean zero and 
variance one and A(/3) := log Eg e~^^ < oo. 

The goal is to study the behaviour of the model as /3 
and d vary. At /3 = the polymer measure is the sim- 
ple random walk, hence the walk is entropy dominated 



and exhibits diffusive behaviour. For /? large the poly- 
mer measure concentrates on paths with low energy and 
the diffusive behaviour is no longer guaranteed. A more 
precisely defined separation between these two regimes is 
given in terms of the quenched and annealed free ener- 
gies: 

FM = lim < lin, l^^i^oW) ^ ,(^). 

Weak disorder corresponds to equality and manifests it- 
self through diffusive behavior. In the opposite, strong 
disorder case the system exhibits localization. The tran- 
sition occurs at a critical /3c- For d = 1 and 2, /3c = 0, 
while < /3c < oo for d > 3 [8] . We only consider d — 1. 

The behaviour of the polymer is described in terms 
of fiuctuation exponents for the polymer endpoint and 
the free energy. For 1 + 1-dimensional polymers there 
are long-standing predictions: For the polymer endpoint 



s{n) 



i^'^ with C = 2/3 for all /3; for the free energy 



fiuctuation exponent VarQ(log Z^(/3)) ~ n?'^, with x — 
1/3. Observe that C, = 2/3 and x = 1/3 satisfy the 
relation x = 2^ — 1, which is expected to hold whenever 
there is localization. 

For the free energy the limiting distribution of the fluc- 
tuations is the Tracy- Widom distribution ^ correspond- 
ing to the Gaussian Unitary Ensemble of random ma- 
trices, that is c{p)[\ogZ^{p) - Fq{P)n)/n^l^ ^ Fgue- 
In fact, a stronger statement holds for the point-to-point 
partition functions. After proper normalization the col- 
lection of partition functions parameterized by the renor- 
malized target point x converges to a universal station- 
ary process called Airy2 10! which has Fgue as its one- 
dimensional distributions. Precisely 



logZ;^(/3;n^/35)- j^^(/3)n 
c(/3)ni/3 



+ a(/?)s2^ Airy2(s), (1) 



where a(/3) is a positive constant which compensates for a 
convex dependence of the free energy on the target point, 
and Z!^(j3; x) is the point-to-point partition function with 
endpoint at x. 
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We introduce the localization length exponent v. 



Varp^_^(s(n)) = £;p^_^ 



{s{n) - £:p^_^[s(n) 



In the weak disorder (diffusive) case v = \/2. In the 
strong disorder case u — Q for typical environments, i.e. 
the polymer is localized. This does not contradict the 
known fact that the variance averaged over environments 
is of order n [llj . For very few environments of total 
probability on the order of n"^/'^, Varp^^ ^ is on the 

order of n"*/^. This explains the apparent discrepancy. 



EXPONENTS FOR INTERMEDIATE DISORDER 

The central focus of this paper is to describe these ex- 
ponents and the limiting distributions in what we have 
termed an intermediate disorder regime, accessed by scal- 
ing /? as we increase the polymer length, 



Pn ■■= Pn 



a > 0. 



(2) 



The a — case is the much studied polymer measure, 
with predictions given on the previous page. For different 
values of a > the "a-polymers" behave in ways that are 
quantifiably different from each other. In this sense the 
intermediate disorder regime is actually an entire family 
of regimes with a as an indexing parameter. The value 
a = 1/4 is critical and is also the most interesting; it is 
described soon. For a > 1/4 the polymer is diffusive and 
behaves as a simple random walk, with -y/n fluctuations of 
the path, i.e. C, = 1/2, v — 1/2. The polymer, rescaled 
by converges to Brownian motion with diffusivity 
constant 1 for almost all environments. Also x = 0. For 
a < 1/4 the exponents C,x, and v depend on a as 

C(«) = ^(1 - = ^(1 - 4a), ^{a) = 2a. 

The leading behaviour of log Z"(/3n^") is of order n^^^". 

The values of C,{a) and x(a) still satisfy x(a) = 
2C,{a) — 1. They linearly interpolate between the KPZ 
[HI III] scalings of C = 2/3, x = 1/3, = at a = 
and the diffusive exponents C = 1/2, x 0,i^ = 1/2 at 
a = 1/4. We now show how the interpolation can be de- 
rived from the values at a = together with the Airy pro- 
cess asymptotics of the point-to-point partition function. 
Let En = Ep^ ^ [s('^)] Slid assume for the moment that 
•n?'^ is of smaller order than En- This assumption means 
that most of the contribution to the fluctuations of the 
free energy comes from paths with an endpoint in a neigh- 
bourhood of En- The logarithm of the number of such 
paths is —E'^jn (by the normal approximation to simple 
random walk) and from (llJ each path contributes en- 
ergy on the order of n^/^^"[Airy2(-E„/n^/^) — Airy2(0)] 



to log Z|^(/3; i?„), so that the log of the point-to-point 
partition function at En is 



n 



,1/3- 



Airy2(-B„/7i' 



2/3N _ 



Airy2(0) 



(3) 



The polymer endpoint favors the value of En which max- 
imizes ([3|. We obtain 



Et 



n 



4/3-a 



Airy2(i;„/n2/3)-Airy2(0) 



It follows that En is of a smaller order than n^/^, and 
since on small distances the Airy2 process is similar to 
the Wiener process we have Airy2(£^„/n^/'^) — Airy2(0) ~ 
^/\En\/ny^- Hence 

El^n'-'^^l 

which gives \En\ ^ Substituting this back 

into Q implies that the fluctuations of the free energy 
are on the order of n^/^('^~'^'^\ which gives the formula 
for x(a). The localization length exponent is determined 
by the condition that the difference between logarithms 
of the point-to-point partition function with endpoints 
separated by the localization length should be of order 
1. For distances \xi — X2 \ ^ n'^^"^ this difference is of the 
order of 



,1/3-a 



Airy2(a;i/n2/3) _ Airy2{x2/n^^^) 



which is of order ^Jx\— X2n^°'- This is of order 1 only 
if i^la) = 2a- We immediately see that for a < 1/4 the 
polymer is localized {v < C)- 

Using the Wiener process approximation the limiting 
distribution for the pair 



'logZ-(/?„;a;n«(")) -logZ,^(/3„;0) ^p-.J^W 



c(q;, /3)nx(") 



C(a,/3)nC(° 



(4) 



is the joint distribution of the maximum value M of 
W{t) — and the point Im where the maximum is 
achieved. Here W{t) is a 2-sided Brownian motion. This 
probability distribution was calculated exactly in |14) . 
The joint density of (tM,M) at {t,a) is g{\t\)ha{\t\)%l)a{,0), 
where g has Fourier transform g{X) = J e''*'^g(s) ds — 
f Ai(z^A)^i, ha has Laplace transform 



/la (A) 



e-^'hais) ds = Ai(pa + ^K0~\ 



and ipai^) has Fourier transform 

V3^(A) = TT^ (Ai(i^A)Bi(iCA + pa) - Bi{i^X)Ai{i^X + pa)) 

with ^ = 2~^^^,p — 2^/3. Here Ai and Bi are the Airy 
functions [15] . 
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C(«) 



11.05 0,1 0.15 0.2 0.25 




For a — 1/4: the localization length is on the same 
order as the displacement of its endpoint (v = ^). The 
polymer is not localized anymore and the latter argument 
breaks down. However the exponents are still correct. 
The scaling behavior in this critical case will be discussed 
in the next section. 

We verify the formulas for C(q^) and x{a) by numerical 
simulations. We fix /3 = 1 and let a vary from to 
.275 in increments of .025. Exponents are computed for 
500 independent Gaussian environments and polymers of 
length n — 5000. See the figures for results. 

CRITICAL REGIME 

In the a — 1/4: case the limiting distribution of the 
polymer endpoint is more involved but still computable. 
From = we do not require any normalization. 

Consider the modified partition function 



Znil^nyX) = Ep 



Y[{1 + I3n~-^u{i, s,))6{s(n) ^ x) 



Expanding the product and then applying the expecta- 
tion with respect to the random walk yields Z„(/3„; a;) — 
ELo/^'^-'/'-^feW, where 

fc 

Here p{i,y) = P(s(i) — y) for random walks starting at 
zero, p'^(ik,Xk) — p{n — ik,x — Xk), and the sum is over 
ordered 1 < ii < . . . < ifc < n with = 0, and free Xj, 
1 < j < fcj with xq = 0. We compute that for fixed k 



n-^I^Jl{x^Jn)/v{n,x^i) ^ 2^1k{x)/ q{1,x), 



where 



lk{x) = \ Q^{tk,X]^)^^W{dti,dXi)Q{ti-ti-X.,Xi~ X.,-x)- 

Here W(i, x) is a Gaussian white noise on x K 
with ¥^[W{t,x)W{s,y)] = S{t - s)5{x - y), g{x,t) = 



exp{-a;^/2t}/v^27rt, gx{tk,Xk) = g{l -tk,x- Xk), and 
the integration is over {0 — to < ti < . . .tk < 1} x M*"' 
with xo — 0. The rT^^'^ term is necessary to keep the 
variance of each term bounded, and the 2^ terms come 
from the local limit theorem for simple random walk. In 
the case k — 1 the convergence is easily seen by a Fourier 
transform computation. We have A(i/3) — — |/3^(l-|-o(l)) 
as /3 — ^ from the mean zero, variance one assumption 
on bj, so 



logE,; 



exp 



^itn i J^{x\/n) /p{n,x\/n)^ 



= -^(1 +o(i))n-^/^EE'^«(^'y)' ^ 

i=l y ^ 

Hereg^(z,y) — p{i,y)p{n—i,x^^—y) /p{n,Xy/n). There- 
fore the limiting distribution is normal with mean zero 
and variance 4(P/^/ti, and it is easily checked that 
2Ii{x) / q{1^x) has the same distribution. For the fc > 1 
terms the convergence is handled by the general theory 
of U-statistics [TB] . 

From these computations we conclude that 
Z^{l3n]Xy/n)/p{n,x^/n) converges to the process 
Zi3{x) := J2k>o('^l^)''^k{x)/g{l,x). The same conver- 
gence also holds replacing Z by Z. Define the process 
A^{x) = logZi3{x) — log g{l, x). It is very important 
that this process is universal, namely it appears as a 
limit for polymer models in the critical regime with 
an extremely general distribution for the random po- 
tentials uj. It is also easy to see that the statistics of 
the density of the scaled polymer endpoint s{n)/^Jn 
is the same as the statistics of the random density 



C/3exp{A^(a 



72} dx. 



Another important property oi Ap{x) is that it interpo- 
lates between a Gaussian process (/3 = 0) and the Airy2 
process (/3 — oo). To see this crossover one has to prop- 
erly rescale the process with /3 to normalize the variance 
of its one-point distribution and its two-point correlation 
function at a fixed distance. The interpolation property 
follows from [17] and [HI [19] , where an exact formula for 
the distribution of Ap(x) was calculated (see also [20j). 
Using the Tracy- Widom formula for asymmetric simple 
exclusion [HI [55] , they obtain: 

^{Ap{x) >s}^j e"''"7(s + 2/?/3-log\/32^-r) dr 
where, 

/(r) = det(/ - K,^)tr ((/ - K,^)-^VMry) 
with K = 2/3^/^, PAiry(x,y) = Ai(x) Ai(y), and 



K„{x, y) ^ I cr(t) Ai(x -I- t) Ai{y + t) dt 
-ir< fx + y" 
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Here a{t) = (1 - e"''*)-! - (nt)-'^ and 

2 /""^ sin (x77 + 77^/12 -g2/,y + 7r/4) 



r3/2 



\A7 



dry. 



In particular it follows that V{Ap{x) < 2P*/^s} 
Fgue{s). It also follows that ^^(a;) is stationary in x 
as its distribution is independent of x. Finally, one can 
show that the process j4^(a;) is continuous and locally 
Brownian, i.e. Aj3{x + 6) — ^/3(a;) is of order ^/S. 

EXACT CALCULATION OF EXPONENTS 

The limiting distribution is also characterized [T71 |TS] 
by observing that is the Wick exponential of a 
stochastic integral. Let B be a one-dimensional Brown- 
ian motion that is independent of the white noise W(t, x). 
Let Eq.x denote expectation over Brownian paths started 
at zero and ending at x at time 1. Then 

»i 



Z^{x) = E, 



'0,x 



:exp:<^ 2^ / W{s,B{s))ds 



The exact density of the polymer endpoint Zp{x)g{l, x) 
in the critical regime is part of the critical polymer mea- 
sure on paths. Define the continuum partition function 



Z^(/3)=Eo 



:exp:<^ 2/3 / W{s,B{s))ds 



(5) 



The polymer measure PJ^ in the critical a = 1/4 case 
scales to a probability measure on paths of length T given 

by 



dP^^(B) =:exp: 



dF{B) 



where P stands for the usual Wiener measure. Write 
for expectation. ^ provides another derivation of 
the fluctuation exponents ({a) and x(q;). For /3 > let 

B{s) := P^B{(3-'^s), W{t,x) := r^W{r'^t, p-^x). 

Scaling properties imply that i? is a standard Brown- 
ian motion and is a standard space-time white noise, 
so via the substitution u = /3*s we have the relation 

P W{s, B{s)) ds = J^^'^ W{u, du and so Z^j, = 

Zfg^j.. From Ya.Tw{ZWrj.) - T^/s we get 



YSiTw{Z^-a 



yf(l-4a) 



which agrees with xi^^) 
E^^ [B{Tf] ^ T4/3 as T 

fixed 13. Since B{T)'^ 



For C(q^) we have 
00 for almost all W and 



/3-2B(/34t)2 we have 



Therefore 



CONCLUSION 

We have identified a new disorder regime for directed 
polymers in = 1 + 1 and computed the corresponding 
wandering exponent C and free energy fluctuation expo- 
nent in addition to a localization length exponent v. 
These exponents vary linearly with the parameter a be- 
tween the KPZ scaling exponents at a = and the simple 
random walk exponents at a = 1/4. At a = 1/4 there 
is a phase transition from intermediate disorder to weak 
disorder. For a > we are able to identify the joint 
limiting distribution of the free energy and the polymer 
endpoint, and give an explicit asymptotic formula for the 
distribution of the point-to-point free energy at a = 1/4. 

The research of all authors is supported by the Natural 
Sciences and Engineering Research Council of Canada. 
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